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As an application of a recently developed variational per- 
turbation theory we find the first 22 terms of the convergent 
strong-coupling series expansion for the ground state energy 
of the quartic anharmonic oscillator. 



Variational perturbation theory yields uniformly and ex- 
ponentially fast converging expansions for any quantum 
mechanical system [1-3]. In this note we adapt this 
theory the problem of calculating strong-coupling expan- 
sions. As an example, we consider the quantum mecha- 
nical anharmonic oscillator and find the expansion co- 
efficients with a high degree of accuracy. The input is 
provided by the exact Rayleigh-Schrodinger perturbation 
coefficients of the ground-state energy up to order 251, 
obtained from the recursion relations of Bender and Wu 
[4]. The resulting strong-coupling coefficients will be cal- 
culated up to the 22nd order with a precision of about 
20 digits. 

The potential of the anharmonic oscillator is 
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V(x) = — x 2 + jx 



(u\g >0). 



(1) 



The Rayleigh-Schrodinger perturbation theory yields a 



power series expansion 



where ef w are rational numbers 
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The superscript BW refers to Bender and Wu. Using the 
symbolic algebra program MAPLE we generate the first 
251 terms exactly. 

The series (2) cannot be used for an evaluation of the 
energy — it has a zero radius of convergence due to the 
factorial growth of the coefficients e^ w . Only for small 
couplings g < 0.1, it yields reasonable approximations if 
truncated at a finite order N, optimally at the integer 
closest to 3/4g . 

The recently developed variational perturbation theory 
[2,3] converts the divergent series (2) into an exponenti- 
ally fast convergent one. The procedure goes as follows 



(see Section 5.13 of Refs. [2,3]). First, the harmonic term 
of the potential is split into a new harmonic term with a 
trial frequency £1, and a remainder: 



SI 2 
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After rewriting 



V{x) = —x 2 + V mt (x), 



with an interaction 
9 



V 1Tlt (x) = ~(rx +x), r = -{l 



9 



SI 2 ), 
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we perform a perturbation expansion in powers of g at a 
fixed r, 



EN{g,r) 
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The calculation of the new series up to a specific order k 
requires only littel additional work, being easily obtained 
from the ordinary perturbation series (2) by replacing u> 
by ^/Sl 2 + gr/2, and by reexpanding (2) in powers of g up 
to the Nth order. This yields the reexpansion coefficients 
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The truncated power series 
W N (g,Q) 



E N [g,-(uj 2 
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is certainly independent of SI in the limit k —> oo. At any 
finite order, however, it does depend on SI, the approxi- 
mation having its fastest speed of convergence where it 
depends least on SI. If we denote the order-dependent 
optimal value of SI by Sljv, the quantity Wn{§, &n) is 
the new approximation to E(g) [5]. In terms of the di- 
mensionless constants 



9 = 9/Q 3 , w=w/fi, 
the approximation can be written as 

W N = (g/g) 1/3 w N (g,io 2 ). 



(10) 



(11) 
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From the approximate energies (9) it is easy to derive 
simple formulas for the coefficients of the strong-coupling 
expansion [3]. We simply expand the function W]y(g,d> 2 ) 
in powers of Co 2 = (g /uj 3 )~ 2 l 3 g 2 l 3 and find 



W N = (ff/4) 1 / 3 
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with the coefficients 



± r ,(2n-l)/3jn) 

n 



9"'" '" w { N '(g,0). 
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Here w^(g,0) denotes the n'th derivative of W]y(g,d> 2 ) 
with respect to Qj 2 at Qj 2 = 0. To calculate these deriva- 
tives, we note that Qj 2 enters the reexpansion coefficients 
(8) in the form 



3=0 



l " J !'" l(4(^-l)/fl)'->. (14) 



The quantities u>jy \g , 0) are therefore given by 
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„bw / (1 - 3j)/2 \ 



(^°) = fe) ne 'b=o = E e " 
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II(/-i-m) (-4/ff)'-'. (15) 



Vm = 



The optimal value of Sljy has the ^-dependence (see 
Ref. [3]) 



il 3 N = gcN (l + 6.85/A 2/3 ) 



(16) 



where the coefficient c is the solution of a simple tran- 
scendental equa- 
tion (c = 0.186 047 272 987 397 512 984 554 740 462. eye- 
ball fit to extremal values of SI. The corresponding order- 
dependent values of the dimensionless coupling constant 
g = g/£i% are inserted into (13) and produce the coeffi- 
cients shown in Table I. Our result for «o agrees to all 
23 digits with the most accurate value for «o available in 
the literature [6]: 

a = 0.667 986 259 155 777 108 270 962 016 919 860 

199 430 404 936 984 060 455 976 663 80. (17) 

For «i and 0.2, our results are consistent with but con- 
siderably more accurate than previous results in Ref. [7] 
and Refs. [8,9] [«i = 0.143 668 800, a 2 = -0.008 627 58]. 

As a further check we have evaluated our strong- 
coupling series at g/4 = 0.1, 0.3, 0.5, 1, 2 (setting u> = 1) 
and compared the numbers with the very precise lower 
and upper bounds of Vinette and Cizek [6]. Table II 
shows that the energies are accurate to about 20 digits 



for all couplings g/4 > 1. The accuracy is limited by the 
precision of the a n . 

Note that our strong-coupling expansion gives good 
results down to very small couplings g — even at g/A = 
0.1, the result agreeing to 7 digits with the bounds in 
[6]. If one wants to work at such small couplings, one 
can easily calculate higher coefficientsa n to increase the 
accuracy. 

In Fig. 1 we show the approach of a n to the asymptotic 
value given in Table I by plotting 
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N 



\( a n) 
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(18) 



on a logarithmic scale. The periodic structure in the data 
is caused by an oscillatory approach of (a n )jv — «n to zero 
(see Fig. 2). The envelope of these oscillations follows 
the curve 
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exp 
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(19) 



For the leading term «o of the strong-coupling expan- 
sion, the behavior (16) suggests a coefficient K\ ps 9.7 
[3]. From the eyeball fits shown in Fig. 1 we extract 
ki « 9.42, 9.05, 8.05, and 7.18 for a n with n = 0, 1, 5 and 
10, respectively. These estimates should be taken with 
some care, however, since we do not know how close we 
are to the asymptotic regime for N ps 65 . . .251, where 
the natural expansion parmeter is 1/ N 1 ! 3 ps 0.25 . . .0.16. 
The TV-dependence of the coefficient ao, for instance, is 
just as well fitted by the dashed line with K\ = 9.7. 

We conclude by mentioning that variational perturba- 
tion theory allows also for a calculation of the imagi- 
nary parts of energies on the left-hand cut in the com- 
plex coupling constant plane [10]. The results not only 
improve the semiclassical values [11] in the tunneling re- 
gime, but also determines the imaginary part at large 
negative g [10,12], where the decay proceeds by sliding 
rather than tunneling. 

W.J. thanks the Deutsche Forschungsgemeinschaft for 
a Heisenberg fellowship. 

TABLE I. Strong-coupling perturbation coefficients a„. 



n 


a n 





0.667 986 259 155 777 108 270 96 


1 


0.143 668 783 380 864 910 020 3 


2 


-0.008 627 565 680 802 279 128 


3 


0.000 818 208 905 756 349 543 


4 


-0.000 082 429 217 130 077 221 


5 


0.000 008 069 494 235 040 966 


6 


-0.000 000 727 977 005 945 775 


7 


0.000 000 056 145 997 222 354 


8 


-0.000 000 002 949 562 732 712 


9 


-0.000 000 000 064 215 331 954 


10 


0.000 000 000 048 214 263 787 


11 


-0.000 000 000 008 940 319 867 


12 


0.000 000 000 001 205 637 215 


13 


-0.000 000 000 000 130 347 650 


14 


0.000 000 000 000 010 760 089 
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15 -0.000 000 000 000 000 445 890 1 

16 -0.000 000 000 000 000 058 989 8 

17 0.000 000 000 000 000 019 196 00 

18 -0.000 000 000 000 000 003 288 13 

19 0.000 000 000 000 000 000 429 62 

20 -0.000 000 000 000 000 000 044 438 

21 0.000 000 000 000 000 000 003 230 5 

22 -0.000 000 000 000 000 000 000 031 4 



The solid straight lines are the best eyeball fits to the enve- 
lope of the data. The dashed line in the «o-data has a slope 
equal to the theoretically expected value 9.7 [3] . 
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FIG. 2. Oscillatory behavior around the exponential ap- 
proach to the limiting value of «o 

TABLE II. Ground-state energies from strong-coupling se- 



ries expansion. The lines labeled "lb" and "ub" are the lower 
and upper bounds of Ref. [6]. 



ff/4 


^max 


E 


0.1 


15 


0.559 146 597 503 562 187 




20 


0.559 146 201 201 805 544 6 




22 


0.559 146 344 373 873 126 9 




lb 


0.559 146 327 183 519 576 3 




ub 


0.559 146 327 183 519 576 7 


0.3 


15 


0.637 991 783 178 536 025 3 




20 


0.637 991 783 171 236 149 3 




22 


0.637 991 783 171 280 381 8 




lb 


0.637 991 783 171 278 528 3 




ub 


0.637 991 783 171 278 529 6 


0.5 


15 


0.696 175 820 765 191 516 9 




20 


0.696 175 820 765 145 887 5 




22 


0.696 175 820 765 145 928 8 




lp 


0.696 175 820 765 145 925 1 




up 


0.696 175 820 765 145 928 5 


1.0 


15 


0.803 770 651 234 273 812 047 6 




20 


0.803 770 651 234 273 769 350 9 




22 


0.803 770 651 234 273 769 354 1 




lb 


0.803 770 651 234 273 756 




ub 


0.803 770 651 234 273 786 


2.0 


15 


0.951 568 472 729 500 011 184 213 69 




20 


0.951 568 472 729 500 011 146 930 27 




22 


0.951 568 472 729 500 011 146 930 52 




lb 


0.951 568 472 729 499 9 




ub 


0.951 568 472 729 500 1 
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